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Abstract 

We study bound states below threshold and resonances above threshold in 
the D°K + and Dfrj systems, using a many-coupled-channel model for non- 
exotic meson-meson scattering applied to states with the quantum numbers 
of cs quark-antiquark vector mesons. We fit the ground state at 2.112 GeV, 
whence the lowest resonances in D°K+ come out at 2.61, 2.72, 3.03, and 
3.08 GeV. The resonance at 2.61 GeV acquires a width of about 8 MeV, 
while its partial P-wave cross section is up to six times larger in D s r] than in 
D°K + , provided a mechanism accounting for Okubo-Zweig-Iizuka-forbidden 
decays is included. The latter finding is in agreement with the observations of 
the SELEX collaboration [1] with respect to the recently reported .D^j (2632) 
resonance. Therefore, we conclude that the -D^j(2632) is probably the first 
recurrence of the Z)|(2112) meson. 
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1 Introduction 



The very fortunate discovery [1] of a new resonance in both the D°K + and Dfrj 
channels has been awaited for over two decades [2]. Excitations of J p = 1~ cs states 
have been foreseen in the past, but with higher masses. For instance, in Ref. [3] 
two states were forecast at 2.73 and 2.90 GeV, one for each of the possible angular 
configurations S and D, whereas in Ref. [4] these states were predicted at 2.773 and 
2.813 GeV. Nowadays, a mass of 2.632 GeV is of no surprise. Just consider the mass 
differences between the ground states and first radial excitations of other J p = 1~ 
mesons, which vary from 0.52 GeV for K* to 0.59 GeV for cc. The 1)^(2632) 
resonance, being 0.52 GeV heavier than the D* ground state, could thus perfectly 
well turn out be the first radial excitation of the D*(2112) meson. However, the 
branching ratio of its two decay modes, i.e., D°K + to Dfrj, with the latter one 
dominating by a factor of six, seems to be problematic [1,5]. Also its decay width 
(smaller than 17 MeV [1]) is smaller than expected for a cs resonance with sufficient 
phase space in a two-meson channel observing the Okubo-Zweig-Iizuka (OZI) [6] 
rule. 

This small decay width, and the puzzling branching ratio to D°K + and D^r], are 
explained by Maiani et al. [7] assuming a dominantly [cd] ds four-quark configu- 
ration, whereas Liu et al. [8] argue that the Dfj(2632) resonance may be a member 
of a flavor 15-plet, which is a mixture of [cu] [us], [cd] ds , and [cs] [ss]. Other 

tetraquark proposals can be found in Refs. [9] and [10], namely [cs] uu+dd and 

[cd] ds , respectively. A radial excitation of the D* is suggested in Refs. [5, 11], 
among other possibilities. 

Here, we shall show that the observed properties are compatible with the predic- 
tions of a multichannel quark-meson model, designed to describe mesonic resonances 
in terms of systems of coupled qq and two-meson states [12]. However, the study of 
non-ground-state D* vector mesons requires a more elaborate approach than the one 
previously applied by us to scalar mesons [13], as we shall explain in the following. 

In the first place, J p = Y~ systems consisting of a charm quark and a strange 
antiquark may appear in two different orbital excitations, namely £ = and £ = 2. 
This implies that instead of one cs channel, as in the case of the _D* (2317) resonance 
[14], we must now consider two such channels. 

Secondly, higher radial and angular excitations of cs can decay into a variety 
of different OZI-allowed two-meson channels, while the influence of nearby closed 
two-meson channels must be taken into account, too. In Table 1, we summarize the 
twenty-eight mesonic decay channels of the cs system that we consider of importance 
in the energy region of interest. 

Thirdly, the phenomenological Resonance-Spectrum- Expansion (RSE) formalism 
employed in Refs. [13, 14] is not practicable in the present case, exactly because of 
the proliferation of relevant two-meson channels involved. In the RSE method, one 
term in the expansion accounts for one radial or angular qq channel coupled to one 
meson-meson channel. Now, each term comes with a different parameter for each 
combination of one of the cs channels and one of the two-meson channels. Thus 
imagining four terms in the expansion, then with two channels in the cs sector 
and twenty-eight channels in the two-meson sector we would end up with many 
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parameters: clearly enough to fit anything, but not very useful for predictions. 

Nonetheless, with a definite choice for the confinement potential, these ambi- 
guities are readily removed. Here, we opt for a harmonic oscillator, which has the 
additional advantage that the relative decay couplings, given in Table 1, can be 
determined through the formalism developed in Ref. [15]. 

Since, moreover, we consider only one flavor state with reasonably heavy decay 
products, it is justified to work in the spherical-delta-shell approximation for the 
3 Pq communication between the cs and two-meson sectors. Thus, we shall employ 
here the model developed in Ref. [16] for cc and bb states, and generalized to other 
flavors in Ref. [17]. The relative couplings are determined with the formalism from 
Ref. [15], where it is assumed that each decay product takes over one of the original 
constituent quarks. Other processes are here classified as OZI forbidden. 



3 



channel 


threshold 




1 J. ' 

relative 


couplings 




GeV 


(L,S) 


to t = 


to I — 2 


D° - K + 


2.358 


(1,0) 


1/72 


1/216 


D + - K° 


2.367 


(1,0) 


1/72 


1/216 


D*° - K + 


2.500 


(1,1) 


1/36 


1/432 


D* + - K° 


2.508 


(1,1) 


1/36 


1/432 


Dt-V 


2.516 


(1,0) 


1/108 


1/288 


d*+ - v 


2.660 


(1,1) 


1/54 


1/648 


D° - K* + 


2.756 


(1,1) 


1/36 


1/432 


D+ - K*° 


2.766 


(1,1) 


1/36 


1/432 


K*(800) - D* 


2.838 


(0,1) 


1/24 


1/216 






(2,1) 




1/216 


D*° - K* + 


2.898 


(1 0) 


1/216 


1/648 






(1,2) 


5/54 


1/3240 






(3,2) 




7/120 


D* + - K*° 


2.906 


(1,0) 


1/216 


1/648 






(1,2) 


5/54 


1/3240 






(3,2) 




7/120 


Df - i 


2.926 


(1,0) 


1/216 


1/648 


D s -<j> 


2.988 


(1,1) 


1/36 


1/432 


d*+ - n' 


3.070 


(1,1) 


1/108 


1/1296 


/ (980) - D* 


3.092 


(0,1) 


1/48 


1/432 






(2,1) 




1/432 


D* s - <P 


3.132 


(1,0) 


1/216 


1/648 






(1,2) 


5/54 


1/3240 






(3,2) 




7/120 


DS(2308) - K* 


3.201 


(0,1) 


1/24 


1/216 






(2,1) 




1/216 


D s * (2317) - 


3.337 


(0,1) 


1/48 


1/432 






(2,1) 




1/432 



Table 1: The pseudoscalar-pseudoscalar, pseudoscalar- vector, vector- vector, and 
scalar-vector real and virtual OZI-allowed two-meson channels, for L = Q 1, 2, 3, 
considered throughout this work, their thresholds, and relative couplings squared to 
each of the two cs channels. For r] and rj' we assume flavor octet and singlet, respec- 
tively. Experimental masses are taken from Ref. [18], except for the _Dq(2308) [19] 
(see also Ref. [14]), and the Kq(S00) (or k) [18]. For the latter resonance, we choose 
the peak mass of 830 MeV from Ref. [13]. 
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2 Inelastic meson-meson scattering 



The 28 x 28 scattering matrix S, as a function of the total center-of-mass energy 
y/s, has the familiar form 



S(y/s) = A; 1 / 2 [1-iK]- 1 [1 + iK] k- 1 ' 2 , 



(1) 



where k and L can be represented by diagonal matrices containing the linear and 
angular momenta of each of the two-meson channels, respectively. The elements of 
the matrix k are determined through the kinematically relativistic expression 



Ask 2 



s - (m u + m 2 i) 2 s - (m u - m 2 



(2) 



where mu and m 2i stand for the meson masses involved in the i-th two-meson 
channel. The 28 x 28 inverse cotangent matrix K, from which scattering phase 
shifts can be obtained, is defined by 



K (y/s) = - [l - k~ l M JXN] 1 [k^M JXJ 



(3) 



In this equation, the diagonal matrices J and N contain the two linearly inde- 
pendent Bessel and Neumann scattering solutions above threshold, respectively, or 
their analytic continuations below threshold, for each of the two-meson channels. 
The diagonal elements of the reduced-mass matrix M are calculated by using the 
relativistic formula 



4s 3/2 Mj (y/s) = s 2 - (m u + m 2i f (mu - m 2i y 



The 28 x 28 matrix X introduced in Eq. (3) is defined as 



*(V5) 



4A 2 V T AV 



(4) 



(5) 



where A is the overall coupling constant related to quark-pair creation, and the 2 x 28 
matrix V contains the relative coupling constants given in Table 1. The diagonal 
2x2 matrix A contains products of the two linearly independent solutions of the 
harmonic oscillator. It is this matrix for which an RSE expansion can be formulated. 

Except for the matrix V, all matrices are diagonal, because we do not consider 
direct interactions in either of the two sectors. More details can be found in Ref. [17]. 

From Ref. [17], we adopt here the effective constituent quark masses m c — 1.562 
GeV and m s = 0.508 GeV, and the oscillator frequency w = 0.19 GeV. However, we 
cannot use the scaled dimensionless radius po from Ref. [17], since it corresponds 
to the maxima of different transition potentials [15], not their effective radii. In 
Ref. [13] we used 5 GeV -1 for in a P wave, which gives po ~ 1) an d 3.2 GeV -1 
for Kn in an S wave, yielding p = 0.66. For simplicity, we choose here po = 0.8, 
as a kind of average value for the different waves under consideration (see Table 1), 
corresponding to a distance of 0.58 fm. The overall coupling A we vary such that the 
ground state of the spectrum, i.e., the D*(2112), gets its experimental [18] mass. 
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3 OZI-forbidden decays 



Our initial findings, with OZI-allowed channels only, are the following. We obtain a 
narrow resonance at about 2.6 GeV. Its width and position depend on the invariant 
radius po. For the above value of po = 0.8, the resonance comes out at 2.613 GeV, 
and has a width of 6.5 MeV. Hence, contrary to what might be expected from naive 
perturbative calculations, we obtain a narrow resonance for a system having more 
than enough phase space for OZI-allowed decay. The reason appears to be twofold. 
First of all, most of the probability of 3 P pair creation goes into closed channels, 
which only give rise to a real mass shift. Secondly, the first radial excitation we are 
assuming for the _Dj}(2632), albeit mixed with the l 3 Di state, has one node (see 
also Ref. [11]), which cannot lie far away from the the delta-shell radius po, ° n the 
basis of general size arguments. In the more realistic picture of Ref. [17], where a 
smooth transition potential was used, the convolution of this potential with the sign- 
changing wave function of the "bare" _Dj}(2632) will lead to a partial cancellation. 
So the total width due to OZI-allowed decays could be considerably smaller than 
what is generally assumed. 

Now, once we have a full solution of the coupled-channel problem, i.e., the com- 
plete S matrix, we may determine branching ratios for other types of decay modes, 
namely OZI-forbidden processes. However, one of these happens to contribute to 
the very Dfrj mode, already included in the list of OZI-allowed decay channels. 
Notice that an OZI decay to Dfi] can only proceed via ss pair creation, whereas the 
non-OZI decay may take place through both uu/dd and ss creation. Precise values 
for such OZI-forbidden decays have to be estimated, since we do not dispose of a lot 
of experimental data on such processes when alternative OZI-allowed modes exist, 
too. Here, we estimate that up to 4-5 percent of the total probability of strong decay 
contributes to OZI-forbidden decays, which we account for in an effective way by 
correspondingly increasing the coupling constant of the OZI-allowed Dfrj channel. 
Note that this estimate of 4-5% seems quite reasonable on the basis of experimen- 
tally known [18] OZI-forbidden two-particle decays of mesons, keeping also in mind 
that the mentioned suppression due to the nodal structure of the wave function may 
not apply to non-OZI decays. In Table 2, we present the results, after refitting the 
overall coupling constant to the D*(2112) mass, for the effect of the non-OZI D^t] 
decay mode on the D*' — > D°K + / D*' — > D s rj branching ratio, as well as the total 
width. We see that the surprising branching ratio D*' -> D s rj / D*' -> D°K + « 6 
observed by the SELEX collaboration [1] is reproduced for a non-OZI contamina- 
tion of about 4.3%, implying a total width of some 8 MeV, which is also compatible 
with the data. Of course, the latter experimental findings need further confirmation. 
What our results show is that such a scenario is plausible. 
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non-OZI (%) 


peak (GeV) 


D*' -> D°K+ 

d*; - d sV 


T (MeV) 


A 





2.613 


10.46 


6.5 


2.67 


0.31 


2.613 


4.47 


6.6 


2.66 


0.62 


2.613 


2.52 


6.6 


2.65 


0.93 


2.612 


1.69 


6.9 


2.63 


1.23 


2.612 


1.18 


7.0 


2.62 


1.85 


2.611 


0.67 


7.2 


2.60 


2.47 


2.611 


0.43 


7.4 


2.58 


4.32 


2.610 


0.17 


7.9 


2.51 



Table 2: The D*' central resonance position, the branching ratio of D*' -^-D°K + to 
D* — > D s i], the total width of the D* , and the effective overall coupling constant, 
as a function of the percentage of non-OZI w.r.t. OZI hadronic decay modes. 
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4 P-wave cross sections for cs decay. 



In this section, we present some figures for the computed partial cross sections, 
corresponding to the case which reproduces the SELEX decay-rate ratio (4.32% in 
Table 2). 

4.1 The D* resonance. 

In Fig. 1, we depict the theoretical cross sections for D°K + and Dfrj in the total- 
invariant-mass (a/s) region of 2.6 GeV. We observe a resonance with a central mass 
about 20 MeV below the SELEX value, and a width of 7.9 MeV. Moreover, we see 




2.60 2.61 2.62 2.60 2.61 2.62 

Vi(GeV) V«(GeV) 

(a) (b) 



Figure 1: Theoretical cross section a for D°K + (a) and D s rj (b) P-wave scattering, 
from a multi-coupled-channel calculation (model parameters taken from Ref. [17], 
except for the transition-potential parameters). 

that about 6 times more D^rj pairs are produced than D°K + pairs. This resonance 
has here vector- meson quantum numbers, i.e., J p = l~. 
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4.2 Further P-wave D°K + cross sections. 

In Figs. 2 and 3 we show other structures we find in D°K + P-wave scattering, in 
the mass region from threshold up to 3.1 GeV. 

Close to the D°K + threshold and well below the -D+77 threshold, we find a 
structure with a width of some 100 MeV (see Fig. 2). Notice, however, that the peak 
cross section is only 0.005 GeV™ 2 . Whether such a weak signal can be observed in 
experiment looks doubtful. This threshold effect is clearly due to the bound-state 
pole of the P*(2112), and does not correspond to a resonance, since we do not 
encounter any additinal pole in this energy region. 




2.40 2.45 2.50 

•v/i(GeV) 



Figure 2: A broad but tiny theoretical signal in D°K + P-wave scattering below the 
D s rj threshold. 
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For energies higher than the D* mass, we find three signals (see Fig. 3): two 
narrow resonances at 2.72 GeV and 3.08 GeV, and one broad signal at 3.03 GeV. 
OZI-forbidden decay modes may very well broaden the two narrow structures. So 




D°K+ 



0.6 



0.4 



0.2 



2.720 2.725 2.730 

Vi(GeV) 




3.00 3.05 3.10 

V^(GeV) 

(b) 



Figure 3: Details of the predicted resonances with central masses at 2.72 GeV (a) 
3.03 and 3.08 GeV (b). Their widths are about 2, 60, and 8 MeV, respectively. 



what exactly can be found in experiment, is not entirely clear from our predictions. 
We only assure that narrow companions to the .0^(2632) should exist. Their po- 
sitions are predicted here within the model error of some 50 MeV. The broader 
resonance at 3.03 GeV, with a width of 50-60 MeV, should be reasonably easy to 
observe experimentally. 
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5 Conclusions 



The experimental confirmation of at least some of the here predicted signals will 
contribute strongly to our understanding of the cs spectrum. 
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